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Abstract. The properties of a system of n = 3 coupled oscillators with linear terms in 
the velocities (magnetic terms) depending in two parameters are studied. We proved the 
existence of a bi-Hamiltonian structure arising from a non-symplectic symmetry, as well the 
existence of master symmetries and additional integrals of motion (weak superintegrability) 
for certain particular values of the parameters. 
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1 Introduction 

The main objective of this paper is to present a study of a system that, although arising from a 
quadratic Hamiltonian of simple form, it is however endowed with several interesting properties 
related with the theory of nonstandard symmetries and the existence of cubic integrals of motion. 

Let us consider the Lagrangian of a system of n = 3 coupled oscillators with an additional 
coupling term linear in the velocities (the so-called magnetic terms) 

L = ^{vl + vj + vl)-{Vi + Vo), 

Vi = b{Q32Vi + qi3V2 + qiiVi), (1) 

where h and k are arbitrary constants, qij denote qi — Qj, and we assume that, for easiness of 
notation, the mass m of the particle is set equal to unity. See Refs. [5] and ^Ij for two recent 
papers dealing with similar systems. 

The starting point is to consider Lq = T — Vq, where T denotes the kinetic energy, as the 
basic original system, and the additional term Vi (linear in the velocities) as introducing a 
deformation. We try to formulate the results in explicit dependence of b in order to study the 
changes of the dynamics when the parameter b varies. 

The Lagrangian L has two exact Noether symmetries generated by the vector fields 
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with associated constants of the motion 

l2=Vi+V2+ V3, 

h = q^2Vi + qi3V2 + q2iV3 - b{qji + gfa + ^13) • 
The first constant Ii is the total energy given by 

h = A(L) -L = ^{vl+vl + vl) + h{ql + + qls) . 

where A = Vi{d/dvi) represents the Liouville vector field. 
The Legendre transformation 

Pi=vi-bq32, P2 = V2-bqi3, P3 = V3 - bq2i 
leads to the following Hamiltonian 

[{Pl + bq32f + {P2 + bquf + {p3 + bq2i)^] + \k{ql^ + qj^ + qh) 

= \{p\ + pI+ pI) + b{Li + L2 + L3) + l{b^ + k) + qh + qh) (2) 

and the above three constants become 

Jl = H, J2=Pl+P2+P3, J3 = Li + L2+L3, 

where Li, i = 1,2,3, denote the components of the angular momentum Lj = qjPk — qkPj- This 
three integrals satisfy { J,., Js} = 0, r, s = 1, 2, 3, and make of H an integrable system for all the 
values of b and k. Thus, we arrive to the following result 

Proposition 1. If we consider the magnetic term Vi as a continuous deformation of the ori- 
ginal system Lq = T — Vq (the value of b representing the intensity of the deformation) then 
this deformation preserves both the existence of Noether symmetries and the integrability of the 
original undeformed system. 

Notice that H can also be written as H = Jq + 6J3 with Jq defined as follows 
Jo = \ {pI +pI+pI)+\ (fc' + k) {ql + ql2 + q^) ■ 

2 Symmetries and integrals of motion 

In differential geometric terms, the dynamics of a time-independent Hamiltonian system is de- 
termined by a vector field on the 2n-dimensional cotangent bundle T*Q of a ?i-dimensional 
manifold Q. Cotangent bundles are manifolds endowed, in a natural or canonical way, with 
a symplectic structure ujq that, in coordinates {{qj,pj); j = 1,2, . . . , n}, is given by 

loq = dqj A dpj, loq = —dOo, 6q = Pjdqj 

(we write all the indices as subscripts and we use the summation convention on the repeated 
index). Given a differentiable function F = F{q,p), the vector field Xp defined as the solution 
of the equation 

i{XF)ujQ = dF 

is called the Hamiltonian vector field of the function F. There are two important properties: 

(i) The Hamiltonian vector field of a given function is well defined without ambiguities. This 
uniqueness is a consequence of the symplectic character of the two-form wq- 

(ii) Suppose that we are given a Hamiltonian H = H{q,p). Then the dynamics is given by 
the Hamiltonian vector field Th of the Hamiltonian function. That is, i(r//)a;o = dH. 
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2.1 Bi-Hamiltonian structure 

At this point we recah that a (infinitesimal) dynamical symmetry of a Hamiltonian system 
(T*Q,uJo, H) is a vector field Y such that it satisfies [y, F//] = 0. On the other hand it is 
known that, in some very particular cases, the Hamiltonian systems can admit dynamical but 
non-symplectic symmetries (for a classification of the symmetries in geometric terms see 
and jni)- In this case we have the following property. 

Proposition 2. Suppose there is a vector field Y that is a dynamical symmetry o/Th but does 
not preserve the symplectic two-form 

vCycJo = U!Y 7^ 0. 

Then (i) the dynamical vector fieldT h is bi-Hamiltonian, and (ii) the function Y (H) is the new 
Hamiltonian, and therefore it is a constant of motion. 

Proof. For a proof of this proposition see 0121 and references therein. A similar property 
is studied in |2Ul I22j for the case of Poisson manifolds. ■ 

Let us denote by I3 the Hamiltonian vector field of the function J3 

i{Y3)uJo = dh. (3) 

Then we have the following property: 

Proposition 3. The vector field Y3, defined by ^ as the canonical infinitesimal symmetry 
associated with J3, can be written as a linear combination of two dynamical but non-symplectic 
symmetries ofTn- 

Proof. We can write Y^ as follows 
>3 = Yi - Yl 



3 ~ -^3 ' 

d d d d d d 

dqi dq2 oq^ dpi dp2 dps 



Y3 = 93— + Qi— + 92— +P3Tr- +PiTr- +P2—, 
3 d d d d d 



dqi dq2 dq^ dpi dp2 dps 



^ c/yi uq2 



Then the vector field Y^ (or Yg ) is neither locally-Hamiltonian with respect to 
£yau;Q = 0, 

nor a symmetry of the Hamiltonian function, Y^{H) ^ 0. Nevertheless it satisfies 
[Ys^,Th]=0. 

Therefore Y^ (or Y^) is a dynamical but non-symplectic symmetry of the Hamiltonian sys- 
tem. ■ 

Thus the dynamical vector field Th is a bi-Hamiltonian system with respect to (u>o,uJa) 

i{TH)uJO = dH, i{TH)uJa = dHa, (4) 

with oja and Ha = Ya{H) given by 

uja = dqiA{dp2 + dps) + dq2A{dpi + dps) + dqsA{dpi + dps), 
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Ha — -J2 — H. 

The function Ha = Ya{H) can be considered as a new Hamiltonian for Th- 
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The bi-Hamiltonian structure (F//, wq, Wa) determines a (1, l)-tensor field R by the relation 
LOa{X,Y)=uJo{RX,Y), yX,YGX{T*Q), 

or, equivalently, R = Qq^oQ^ [with Q{X) = i{X)uj G A^(r*Q)]. We have obtained 

d d d 

R = -7^(^{dq2 + dqs) + -7^(S){dqi + dqs) + t— (8)((igi + dq2) 
oqi oq2 oqs 

+ i^^idp2 + dps) + T^^l'ipi + dp3) + T^(8)(a!pi + dp2). (5) 
dpi ap2 dps 

Starting with the basic Hamiltonian system (ujqjTq = Tn^dHQ = dH), and iterating R, we 
can construct a sequence of 2-forms tOk, vector fields F^, and exact 1-forms dHi^., k = 1, 2, . . ., 
defined by cD^ = cDq o i?'^, F^ = i?'^(Fo), and dH^ = R^{dHQ), that leads to sequence of bi- 
Hamiltonian vector fields. The two first steps become 

i{To)uJi = z(Fi)u;o = dHi, i{To)uJ2 = «(Fi)u;i = i(F2)u;o = dH2, 
where u!2, F2, and H2, are given by 

L02 = 2u;o + uJu F2 = 2Fo + Fi, H2 = 2Hq + Hi. 
So, in this particular case, the sequence {uJk,Tk,dHk) closes over itself after two steps. 

2.2 Master symmetries 

The function T = T{q,p) is said to be a generator of constants of motion of degree m if it is not 
preserved by the dynamics but it generates an integral of motion by time derivation 

d r/™~^ J™ 

dt ^ ' ' dt™-! ^ ' dt™ 

In differential geometric terms, a vector field X ot].T*Q that satisfies 

[Fh,X]/0, [Th,X]=Q, X = [Th,XI 

is called a "master symmetry" or a "generator of symmetries" of degree m = 1 for Th (51 151 171 [T7]. 
If a Hamiltonian H admits two of such generators, Ti, T2, (master integrals) then the function 
I12 defined by 

I12 = hTi - I1T2, Ir = ^Tr, r = 1,2, 

dt 

is a new constant of motion (independent of Ji and 12). In geometric terms, if we denote by Xr, 
the Hamiltonian vector fields of T^., then X^ is Hamiltonian as well, and we have 

i{Xr)uJo = dTr, i{Xr)ujo = dfr, = F//(r,,), r = 1,2, 

so that the vector field X12 defined as 

X12 = T1X2 - T2X1 + faXi - fiX2, i(Xi2)c^o = 

is a symmetry of the system [T*Q,ujq,H) with the function I12 as its associated constant of 
motion. 
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The Hamihonian vector field Xi of the function Ti = gi + (72 + 93 satisfies 

i{Xi)u;o = dTi, [Th,Xi]=Xi, = — + — + — , [Th,Xi]=0, 

dqi dq2 dq^ 

so Xi is a "generator of symmetries" of degree m = 1 for Similarly if we consider T2 = 
qiPi + i?2P'2 + 93^3 we obtain 

i{X2)u;o = dT2, [Th, X2] = X2, {rH,X2] = 0, 

where the asterisk means that the Lie bracket vanishes only if A; + 6^ = 0. Hence in the particular 
case /c + 6^ = 0, the system is superintegrable with an additional integral of the form I12 that 
becomes 

J4 = h2 = Li{p2 - P3) + L2{p3 - Pi) + Lsipl - P2). (6) 



2.3 Third-order integrals 

Linear integrals of motion be obtained from exact Noether symmetries and quadratic integrals 
from Hamilton-Jacobi (Schrodinger) separability but third-order integrals must be obtaining by 
other alternative procedures; because of this, the number of known integrable systems admitting 
cubic in the momenta constants is very limited |H1 1101 1131 1211 123| . For n = 2 we can point out 
the Fokas-Lagerstrom and the Holt potentials for the Euclidean plane |S1 112j , and three of the 
Drach potentials for the pseudoeuclidean plane ^1 E| ■ For n = 3 two very special systems are 
well known: the three-particle Toda chain and three-particle Calogero-Moser system which are 
both endowed with a cubic constant of the form 

J = - [pf + P2 + P3) + terms of first order. (7) 

At this point we recall the following property: polynomial constants of the motion for mechanical 
systems (Riemmanian metric with a potential) are either even or odd in the momenta (for 
a discussion of this result see, e.g., ^2)- Therefore if J is a cubic integral for a Hamiltonian 
such as the Calogero-Moser or the Toda systems then it must also contain linear terms but 
neither quadratic nor independent ones. Nevertheless this property is not true for the more 
general case of a non-mechanical Lagrangian. In our case, due to the presence of linear terms in 
the velocities (magnetic terms) when looking for third-order integrals we must assume for the 
function J4 a polynomial expression containing also even powers. 
We start with a third-order function J of a general form 

J= ^ aijkP\P2P?,+ X] bijkPlP2P3+^(^iPi + d, 
i+j+fc=3 i+j+k=2 i 

where aijk, hjk, Ci, and d are functions of qi, q2, qs- Then the condition {H,J} = leads to 
a system of equations; some of these equations restrict the expressions of the coefficients aijk, 
i + j + k = 3, of the higher-order terms; the other equations couple the derivatives of Vi 
and Vo with the coefficients aj^fc, bijk, Ci, and d or with the derivatives of these coefficients. 
Unfortunately, the general case is difficult to be studied but some particular cases can be analyzed 
assuming simple particular forms for the coefficients ttyk- 

We have studied the existence of a cubic integral of such a particular form for this three- 
particle system and we have arrived (we omit the details) at the existence of the following 
function 

-'4 = ^ {pI +pI+ pI) - b{q32pl + qiwl + mpl) 

+ 52 [zipi + Z2P2 + Z3P3) + 9b^ {qli + g?3 + ^32) > (8) 
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where zi, Z2, and z^, are given by 

zi = ^{qli + qls) + qh, Z2 = 4:{ql2 + q2i) + qls, ^^s = 4(^13 + gfa) + 921- 

The important point is that the function J4 is a constant of motion, not for all the values of 
the parameters (6, fc), but only in the particular case of k and b satisfying k = 86^. We have 
verified that in this case the four functions H, J2, J3, J4 are independent and that the system 
(J2, H, J4) is involutive (nevertheless {J3, J4} 7^ 0). Thus, when k and b are related hy k = 86^ 
then the system becomes superintegrable with an additional third-order integral. 

We finally note that J2 is linear, H quadratic, and J4 cubic, and that they have expressions 
of the form 

J2=Pi+P2+P3, H =pI+pI+pI + --- , J4=pI+pI+pI + --- (9) 

that closely resemble the expressions of the three first constants of Toda or the Calogero-Moser 
systems. 

3 Final comments 

The Lagrangian system (Hamiltonian (j2j) is a system endowed with a bi-Hamiltonian struc- 
ture, non-symplectic symmetries and master symmetries. Moreover for a particular value of the 
parameter b, it becomes superintegrable with a cubic integral. 

The non-symplectic symmetry is 6-independent and hence the bi-Hamiltonian structure 
{TH,uJo,tLia) is the one endowed by the the original oscillators system Lq = T — Vq that is 
preserved under the 6-deformation introduced by the magnetic term Vi. Concerning the master 
symmetries, the first one is always present but the second only exists if A; + 6^ = 0. So in 
this particular case the system becomes (weakly) superintegrable. This particular case is very 
peculiar because the Hamiltonian ^ becomes a system of n = 3 particles with velocity-linear 
couplings but without the potential term. 

Finally the particular k = 86^ case is very interesting. First because, as stated above, the 
number of known integrable systems with third-order integrals is very small. Second because 
the existence of this property for a very particular case of the parameters resembles situations 
characterizing other known systems (e.g., the Henon-Heiles system is only integrable for some 
very particular values of a parameter). Finally because the expression obtained for J4 in (jHJ 
shows similarity with the third-order integral of the Toda or the Calogero-Moser systems. 
This leads to the necessity of the study of the more general n > 3 case, that poses the question 
of the nature of the interaction (in the Toda chain the interactions are between neighbours but 
in the Calogero-Moser case they are between every two particles). On the other hand this also 
leads to the search of a Lax pair generating the three functions (J2, H, J4), as traces of the 
powers of an appropriate matrix. Nevertheless, notice that this Lax structure (if it exists) must 
correspond to the L system with the additional condition k = 86^; so it cannot be obtained as a 
continuous deformation of a previous Lax pair for the Lq system. We think that these are open 
questions that must be investigated. 
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